arXiv: 1506.07102vl [quant-ph] 20Jun2015 


Controlling the nonclassical properties of a hybrid Cooper pair box system and an 
intensity dependent nanomechanical resonator. 

C. Valverde,i-2’3>Q V.G. Gongalves,^ and B. Baseia^ 

^ Campus Henrique Santillo, Universidade Estadual de Goids, 

Rod. BR 153, km 98, 75.132-903 Andpolis, GO, Brazil 
‘^Universidade Paulista (UNIP), Rod. BR 153, km 1, 74-845-090 Goidnia, GO, Brazil 
^Eseola Superior de Negoeios-ESUP, Av. Antonio Eidelis 515, 74-840-090 Goidnia, GO, Brazil 
^Instituto de Pisiea, Universidade Pederal de Goids, 74-001-970 Goidnia, GO, Brazil 

(Dated: June 24, 2015) 

We employ a more realistic treatment to investigate the entropy and the excitation-inversion 
of a coupled system that consists of a nanomechanical resonator and a superconducting Cooper 
pair box. The procedure uses the Buck-Sukumar model in the microwave domain, considers the 
nanoresonator with a time dependent frequency and both subsystems in the presence of losses. In¬ 
teresting results were found for the temporal evolutions of the entropy of each subsystem and of 
the excitation-inversion in the Cooper pair box. A comparison was also performed about which of 
these two subsystems is more sensitive to the presence of losses. The results suggest that appro¬ 
priate choices of the involved time dependent parameters allow us to monitor these two features of 
the subsystems and may offer potential applications, e.g., in the generation of nonclassical states, 
quantum communication, quantum lithography. 

PACS numbers: 42.65.Yj; 03.65.Yz; 42.50.Nn 

I. INTRODUCTION 

In the last years the investigations on nanomechanical 
systems [T] have rapidly been developed. The rush in this 
direction was estimnlated by various perspectives and ap¬ 
plications before unsuspected. The enormous progress in 
the research of nanomechanical systems has also placed 
the dream of controlling the interface between the quan¬ 
tum and classic worlds in a realistic way. This was shown 
by the emergence of hybrid quantum systems mm, which 
are intended to achieve a coherent transfer of quantum in¬ 
formation from a single quantum emitter (e.g., supercon¬ 
ducting qubits, cooper pair box, microwave resonators, 
quantum dots, etc) and a solid state mechanical res¬ 
onator, which moves in quantum physics giving birth to 
a new paradigm. In addition, these systems have re¬ 
ceived considerable attention due to their diverse poten¬ 
tial applications, including metrology (mass force or spin 
ultra-sensitive detectors), based on the remarkable sens¬ 
ing properties of nanomechanical resonators and their 
very low mass and lightness. 

An important focus of quantum optics is concerned 
with the atom-photon system. Inspired on the various 
tests applied to this coupled system and on the several re¬ 
sults obtained, including the limitations, the researchers 
have passed from the light domain to the microwave do¬ 
main of the superconducting version, the quantum elec¬ 
trodynamics circuit. This system furnishes a new test for 
microwave '‘’‘photons’^ that interact with superconducting 
qubits [HIS]. In this scenario the atom is substituted by a 
Cooper pair box (CPB) while the photon is substituted 
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by a nanomechanical resonator (NR). So, the atom-field 
interacting system goes to the CPB - NR interacting 
system with the concomitant passage from the optical 
domain to the microwave domain. 

There are few works in the literature that treat the 
interaction between a CPB and a NR when the later 
has either a time dependent frequency [SI] or a time 
dependent amplitude of oscillation [TOl [TT]. The well 
known Jaynes-Cummings model (JCM), which describes 
the interaction of a single two-level atom and a single 
mode of a quantized radiation field, is the simplest model 
for this system and provides exact solutions. Analogous 
to the CPB - NR system, but since the year 1963, 
many others studies have previously implemented the 
JCM to describe the atom-field interaction [HI i9i 
nn. Some generalized models were also constructed and 
extensively studied [TSHIS] . The examples include the 
study of the mentioned systems in presence of the Stark 
effect [T9l . e.g., to investigate quantum nondemolition 
measurements [2TH24]- Usually the investigations assume 
the field initially in a (pure) coherent state. As the atom- 
field interacion is turned on, the field state changes with 
time the evolution. Then, if its coherence is lost, the field 
state becomes nonclassical [25H2I]- Some references to 
this subject are, e.g., on nonclassical properties of a state 
[28H3Q] : generation of superposition states [31]|32]; on the 
degree of nonclassicality of a state [33] and: sculpturing 
coherent states to get Fock states [34] - plus references 
therein. 

The traditional JCM was extended to the case of 
intensity-dependent eoupling, proposed by B. Buck and 
C.V. Sukumar [35] in order to study the infiuence of the 
field intensity, via its excitation number, upon the atom- 
field system, where a single atom interacts with a single 
mode of an optical field. This model was generalized by 
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V. Buzek [36] to include a new coupling, with time de¬ 
pendent intensity. In the present work we will employ 
the model by Buck-Sukumar (BS) to study the CPB - 
NR system, namely: we will suppose the coupling being 
dependent of the intensity of the NR oscillations and also 
that it changes with time, as assumed in [36]. In addition, 
we consider a more realistic scenario with the presence 
of dissipation effects and verify in which way they af¬ 
fect the excited level of the CPB, the excitations of the 
NR, and the dynamical properties of the entire CPB - 
NR system. Some points considered here are: how dis¬ 
sipation spoils the system operation and in which way 
the detuning could prevent it, allowing us the control of 
entanglement features, collapse-revival effects, and oth¬ 
ers. The results obtained indicate the possibility of some 
potential applications [37H42] . 


II. THE HAMILTONIAN SYSTEM 

A superconductor CPB charge qubit is adjusted to 
the input voltage Vi of the system, through a capacitor 
with an input capacitance Ci. Following the configura¬ 
tion shown in Fig. [^we observe three loops: a small loop 
in the left, another in the right, and a great loop in the 
center. The control of the external parameters of the sys¬ 
tem can be implemented via the input voltage Vi and the 
three external fluxes <I>^ and The control of the¬ 
ses parameters allows us to make the coupling between 
the CPB and the NR. We consider H = 1 and assume as 
identical the four Josephson junctions of the circuit sys¬ 
tem, having the same Josephson energy Ej; the external 
fluxes ^L and are also assumed as identical in magni¬ 
tude, although they have opposite signs = —^r = 
(see Ref. i)- So, taking into account the decay in the 
excited level of the CPB and dissipation in the NR, we 
can write the total Hamiltonian of the system as follows. 


H = uj{t)a^a + ^ujc{t)d-z + A(t) (^aV dtd(j+ -h V — ij{t) |e)(e| — i6{t) 


a^d. 


( 1 ) 


In the Eq. 0 the first term describes the NR, the sec¬ 
ond describes the CPB, the third represents the intensity 
dependent interaction introduced by the BS model, here 
including the time dependent parameter A(t) introduced 
by Buzek, the fourth term stands for the time de¬ 
pendent loss affecting the CPB and the fifth term S{t) 
stands for the time dependent loss that affects the NR. 
In the Eq. Q d^ (d) is the creation (anihilation) operator 
for excitations in the NR, (j+ (d_) is the raising (lower¬ 
ing) operator for the CPB, and is the z-component of 
the Pauli spin operator. In the Hamiltonian of Eq. 0 
stands for a non Hermitian Hamiltonian {NHH). Usu¬ 
ally, quantum mechanics works with hermitian Hamilto¬ 
nians; however, one can found a lot of papers in the lit¬ 
erature using NHH : one of them appears in Ref. [44] . 
used by the authors to map a wave guide system, metal- 
silicon, where an optical potential is modulated along 
the length of the wave guide to get a non reciprocal light 
propagation; another application appears in Ref. [45] . 
employed by the authors to show occurrence of entangle¬ 
ment in many body systems; in the Ref. [46] the author 
study optical realizations of a relativistic NHH; in Ref. 
El the authors use a NHH and a convenient algoritm 
to generalize a conventional theory; in Ref. |48| a NHH 
is applied to get information about input-output chan¬ 
nels; an approach in Ref. m neglects the hermiticity 
to study canonical transformations in quantum mechan¬ 
ics; Ref. m applies a NHH to solve a quantum master 
equation; Ref. [51] uses a NHH in a new approach to 
study a weak spectral density of H-bonds and attenua¬ 
tions in parts of the system; Ref. [52], studies a NHH to 
show that it exhibits real eigenvalues; Ref. [53] employs 


a canonical formulation to study a dissipative quantum 
mechanics that exhibits complex eigenvalues; Ref. [54] 
investigates a NHH in non commutative spaces; in the 
Ref. [55] a NHH is considered to include damping in a 
system described by JCM and to study entangled states 
of atoms in distinct cavities; etc. 

In addition to the previous references concerned with 
NHH^ here it is pertinent remembering an alternative 
treatment by Carl M. Bender [56], who defended the effi¬ 
cacy of the non hermitian PT-symmetric Hamiltonians. 
Accordingly, although the hermiticity of the Hamiltonian 
is sufficient to guarantee the two essential properties of 
quantum mechanics (namely, reality of eigenvalues and 
unitary of the time evolution), it is not necessary. Actu¬ 
ally, various complex Hamiltonians are also able to guar¬ 
antee these two properties and they include the real sym¬ 
metric Hamiltonian as a special case, required by the PT- 
symmetric quantum mechanics. (Indeed, even the sup¬ 
posed guarantee assigned to hermitian operators is con¬ 
troversial isaisH].) Although the theoretical approach 
in this line began around 1998, the first experimental 
results came in the last five years [59H64] . 

In the Hamiltonian of Eq. 0 <^+ = |e)(fl|, <3-- = |5)(e| 
and dz = |e)(e| — \g){g\ stand respectively for the transi¬ 
tion ((j+,(j_) and atomic inversion (dz) operators; they 
act in the space of atomic states with frequency Co’c(t); 

and d are the (non hermitian) creation and annihilla- 
tion operators of the NR excitations, with frequency 
00 (t) = ooq E f{t): f{t) is a t-dependent function, which 
can also be constant; is the decay parameter of the 
CPB, from its excited level |e) to its fundamental level 
1^); S{t) is the decay parameter of the NR and A(t) is 
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FIG. 1. Model of a Cooper pair box system coupled to a nanoresonator. 


the coupling coefficient between the CPB and the NR. 

As known in the literature the coupling parameter A(t) 
can be written in the form [65] , 

p' 

where the quantization volume V (t) is dependent on time 
and takes the form V (t) = is the permissiv- 

ity constant; dy = e{i\f\j) is the matrix element of the 
dipole between the two CPB states \i) and \j); e is the 
elementary charge, r is the position vector and, \i) and \j) 
play similar roles of the atomic states |e) and \g) in the 
scenario of atom-field interaction. The notation |e), \g) 
will be used from now on. It stands respectively for the 
excited and ground states of the CPB. 

Substituting the expressions of the field frequency uj {t) 
and quantum volume V (t) in Eq. ([^, we can write it as 


A(t) = Ao 



( 3 ) 


where Aq = The control of the parameters 

uj{t) and A(t) is provided by an external field that acts 
upon the NR and changes the magnetic flux (cf. Fig. 

0 ' 

The wave function that describes the CPB - NR sys¬ 
tem as function of the time t can be written in the form. 


CX3 

1^ (^)) = E W |e, n) + Cg,n (t) \g, n)], (4) 


where the coefficients Ce,n (^) and Cg^n (^) stand respec¬ 
tively for the probability amplitudes to find the entire 
system in the states |e,n) and \g^n). This notation rep¬ 
resents the CPB in its excited state |e) and fundamental 
state 1^) with n excitations in the NR. We will assume 
the subsystems CPB and the NR are decoupled at t = 0 
with the CPB in its excited state |e) and the NR pre¬ 
pared in a coherent state ja), this later expressed by the 
superposition of Fock states. 


|a) = F 1^) ’ 

n=0 


( 5 ) 


where In this way the wave function of 

the whole system in the the initial state can be written 
as |T(0)) = \e)\a) = |e, n), since the initial con¬ 
ditions restrict the probability amplitudes to = 0 

for all values of n = 0,1, 2, 3,... and |C'e,n(0)|^ = 1- 

By analyzing the time evolution of the CPB - NR sys¬ 
tem, described by the time dependent Schrodinger equa¬ 
tion. 


ih 


dt 




(6) 


and the Hamiltonian H given in the Eq.([^, we find the 
set of equations of motion for the coefficients Ce,n (^) and 


dCe,n jt) 
dt 






tW - 


2 


n5{t) Ce,n {t) - iX{t){n + l)C^,n+l {t) 


( 7 ) 















—i{n + l)cj(t) + i —— (n + l)S(t)^ Cg^n+i (t) — iX{t){n + l)Ce,n (t). 
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dCg^n+l {t) _ f 
dt ~ \ 


( 8 ) 


Solving the set of equations Q and ^ we obtain the 
solutions for Ce,n{t) and To this end, we have 

used the Runge-Kutta method of 4^^ order; with theses 
solutions we determine the quantum dynamical proper¬ 
ties of the system, including the entanglement that affects 
the subsystems CPB and NR. 

III. ENTROPY OF THE CPB - NR SYSTEM 

Recently, several authors have employed various meth¬ 
ods to investigate the dynamics of entanglement [HEa- 
[ 68 ] . Here the name “entanglement” means “mixing of 
states”, whose degree can be measured by the entropy. 
The Von Neumann entropy offers a quantitative measure 
of the system disorder or the degree of a quantum state 
purity, e.g., as shown by Phoenix and Knight [69]. Here 
the entropy defined in the form, 

5 ' nr ( cpb ) = —TrN(c){PN(c) ^^Pn(c))^ ( 9 ) 

is a measure related to the entanglement of one of the 
two interacting subsystems, with Pn{C) = Trc{N){PNc) 
where pN{c) is the density operator that describes the 


state of the subsystem NR(CPB) and p^c is the same 
operator for the entire system CPB - NR. The quantum 
dynamics represented by the Hamiltonian in Eq. 0 cre¬ 
ates the entanglement in both subsystems of the CPB 
- NR. In the following discussion we use the quantum 
entropy of Von Neumann as a measure of the degree of 
entanglement. For a quantum system with two compo¬ 
nents the entropy obeys the Araki-Lieb theorem, which 
states that |>Scpb ~ 'S^nrI ^ ^ >5cpb T >5nr- One 

consequence of this inequality is that, if the total system 
is initially prepared in a pure state, then the entropies of 
the components of the system remain with equal value in 
their subsequent time evolution. Then it is sufficient to 
study one of them to know both. 

So, assuming the two subsystems in pure states at t = 0 
the entropies of the CPB and NR subsystems, which 
are found via the BS model, are identical: 5'cpb(^) = 
5'nr(^)- For example, the entropy of the NR is found 
from the equation, 

5 'nr(^) = ~ In (^Nr(^)) + '^Nr(^) ('^Nr(^))] 

( 10 ) 

where : 


^ / oo oo ^ 

^NRit) = o E + E ± ^ 


+4 


n=0 

oo 


n=0 




^ n=0 


n=0 




n=0 


l/ 2 x 


( 11 ) 


The results obtained from the calculations are shown 
in the Figs. (|^|^|4]and|5| 


IV. THE CPB EXCITATION INVERSION 


The CPB excitation inversion, denoted by /(t), is an 
important observable of the two-level systems. It is de¬ 
fined as the difference of probabilities of finding the CPB 
in the excited and fundamental states. The mathemati¬ 
cal expression representing this property is. 


m = J2[\CeAt)f-\C,,n+l{t)f 

n=0 


( 12 ) 


The results achieved from the Eq. ( [T^ , for various values 
of parameters, are shown in Figs. ^^and[^. 


V. RESULTS AND DISCUSSIONS 

Firstly, for f{t) = 0 we will consider the time evolu¬ 
tion of the NR entropy for various values of the decay 
coefficients, 7 (t) and 6{t). We assume the NR initially in 
a coherent state \a) with average number of excitations 
(h) = \a\‘^ = 9 and the resonant case ooq = = 2000Aq; 
other values of parameters are used in Fig. In an ideal 
system the parameters 7 (t) and 6{t) are null, as used in 
Fig. (a). In this figure, the maximum of the NR en¬ 
tropy is close to In 2 . After the start of the interaction the 
NR entropy gradually goes to its minimum, then returns 
to its maximum and remains oscillating regularly. 

For small values of decay in the CPB, as 7 = O.OOIAq, 
and an ideal NR {6 = 0) the maximum value of the 
entropy shows no significative changes for small times 
(cf. Fig. [^(b)). For larger values of time the amplitude 
of the entropy oscillations diminishes while the entropy 
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FIG. 2 . Time evolution of the entropy for different values of 
the parameters 7 (t) and S{t) for: (n) = 9, cjq = c^c = 2000Ao, 
/(t) = 0. (a) 7 = 0,0Ao and 6 = O.OAo; (b) 7 = 0,001Ao and 
(5 = 0.0Ao. 


FIG. 3. Time evolution of the entropy for different values of 
the parameters 7 (t) and S{t) for: (n) = 9, cjq = c^c = 2000 Ao, 
f{t) = 0 . (a) 7 = O.OOlAo and 6 = O.OOlAo; (b) 7 = O.OOlAo 
and 6 = O.OIAq. 


itself moves away from the value S{t) = 0 . If instead we 
include a small decay in the NR, as S = O.OOlAo, the 
entropy changes as follows: although the amplitude of 
the entropy oscillations again decreases, now the entropy 
moves slowly to its minimum value S{t) = 0 (cf. Fig. 
(a)). For larger value of S{t) the entropy movies rapidly 
to zero, which is due to the passage of both subsystems 
to their respective ground states. In this case the entropy 
loses its periodicity (cf. Fig. Kb)). 

Here we consider the nonresonant case, namely, for 
f{t) = A = const^ with A ^ ujc^ ^o- Comparing the Fig. 

(a) with Fig. (b), both have the same values of the 
decay coefficients 7 (t) and 5{t)^ we observe the amplitude 
of entanglement decreasing as the detuning A increases, 
i.e., the detuning turns the entropy oscillation greater 
while destroying its periodicity. In this scenario (A 7^ 0) 
if the CPB decay parameter ^(t) increases the entropy 
goes quickly to zero (not shown in figures). Comparing 
the Fig. [^(b) with Fig. [^(a) we see a similar behavior of 
the average values of the entropies whereas the amplitude 
of oscillations diminishes in the case A 7 ^ 0 . All these 
comparisons show that the infiuence of detuning upon 
the NR entropy is greater in the ideal NR ((5 = 0). 
In the resonant case, when the decay parameters 7 (t) 
and 5{t) increase the NR entropy tends to zero. As a 
consequence of the mentioned Araki-Lieb theorem, the 
same occurs in the CPB. 

Let us now consider the variation in the detuning pa¬ 
rameter A: we fist take A 7 ^ 0 and /(t) = 7 ^ sin(co’/t), 
where 77 and ujf are parameters that modulates the NR 
frequency. Our discussion is limited to the condition 



FIG. 4. Time evolution of the entropy for different values of 
the parameters ^{t) and 5{t) for: (n) — ujq — ujc — 2000Ao 
(a) 7 = O.OOlAo and 6 = O.OAo, fit) = lOAo; (b) 7 = O.OOlAo 
and 5 — O.OOlAo, fit) — 20Ao. 


T] ujc, ojQ and also assuming that uof is small to avoid in¬ 
teraction of the CPB with other modes of the NR. We 
have chosen various values of amplitude modulations 77 
to verify the entanglement properties between the CPB 
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and NR. We use various values of frequency modulation 
uof to see its influence upon the CPB - NR entanglement 
(cf. Figs. [^. 




FIG. 5 . Time evolution of the entropy for different values of 
the parameters 7(t) and S{t) for: (n) = 9 , c^o = c^c = 2000Ao, 
7 = O.OOlAo and 6 = O.OOlAo; (a) 77 = 20Ao, oof = Ao; (b) 
T] = 2 OA 0 , ojf — 2 OA 0 . 


Comparing the Fig. (b) with the Figs. m (a), m 
(b) we note the entropy exhibiting periodic and quasi 
periodic oscillations: when the parameter 77 increases the 
amplitude of these oscillations diminishes. This means 
that the modulation of the NR sinusoidal frequency is 
important to stabilize entanglements in the CPB - NR 
system. 

Next we consider the time evolution of the excita¬ 
tion inversion {El) for different values of the coefficients 
7 (t) and 5{t). In the resonant case, f{t) = 0, we have 
fixed again the NR with average number of excitations 
(n) = |(ap = 9 and ujq = ujc = 2000Aq to calculate the 
El of the CPB (cf. Figs. (a), (b) and (c) ); we 

see in these three cases the El exhibiting similar collapse- 
revival effects, but different amplitudes; these amplitudes 
diminishes for larger values of the parameters 7 (t), S{t). 
In case of detuning, with f{t) = A = const.^ A ^ cJcCJo, 
Fig. [^(a) shows absence of the collapse-revival effect; in¬ 
deed even the El is also absent since these oscillations 
of the NR excitation characterize no inversion, as ex¬ 
plained in the “alert” below. On the other hand, for 
f{t) = r]sm{ujft) the Fig. (b) stands for a variable de¬ 
tuning with maximum value 77 , equal to the value of the 
fixed value 77 = A = 20, Fig. (a). For small times 
the system presents the collapse-revival effect, whereas 
for large times this effect vanishes, with concomitant at¬ 
tenuation in the EL In the (c) we use a large value 
for uf and the same maximum value 77 = 20 employed 




(b) 



JL JL JL J- -L J- -1- /I 




(c) 

FIG. 6 . Time evolution of the excitation inversion for different 
values of the parameters 7(t) and S{t) for: (n) = 9 , Ct;o = 
ujc = 2OOOA0 (a) 7 = O.OOlAo, S = O.OAo and f{t) = 0 ; (b) 
7 = O.OOlAo, S = O.OOlAo and f{t) = 0 ; (c) 7 = O.OOlAo, 
6 = O.OO5A0 and f{t) = 0 . 


in the Fig. (b): here we observe no occurrence of the 
collapse-revival effect, whereas the amplitude of oscilla¬ 
tions decreases for long times which is due to the losses 
affecting both subsystems. This later result differs from 
the case /(t) = A, as shown in Figs. |^(a),[^(b) and|^(c). 
Alert: two different properties must be distinguished in 
the Fig. [^(a): one of them is the excitation inversion, the 
other is the oscillation of excitation. The first is identified 
when the oscillations changes the signs (+) ^(—) suc¬ 
cessively; the second neglects these signs. Ignoring these 
aspects can confuse fluctuations with inversion. For ex¬ 
ample, Fig. (a) shows oscillations, but not inversion; 
then, what is really shown in this figure is a collapse of 
oscillations, not a collapse of the excitation inversion. 
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(b) 



(c) 

FIG. 7. Time evolution of the excitation inversion for different 
values of the parameters 7 (t) and S{t) for: (n) = 9, c^o = 
a;c = 2 OOOA 0 , 7 = O.OOlAo, (5 = O.OOlAo (a) f {t) = 20Ao; (b) 
r] = 20Ao and cj/ = Ao; (c) 77 = 20Ao and ujf = 20Ao; 


VI. CONCLUSION 

We have studied the dynamical properties (entropy 
and excitation inversion) of an interacting system com¬ 
posed by a CPB and a quantized NR. We have assumed 
the CPB initially in its excited state |e) and the NR 
initially prepared in a coherent state \a). We also as¬ 


sumed the whole system CPB - NR described by the 
BS model via single excitation transition. Concerning 
the entropy, which is connected with entanglement (mix¬ 
ture) of states, we have studied its time evolution in the 
presence of loss in both subsystems, for the resonant case 
(A = 0) and nonresonant case. The influence upon the 
entropy of a (sinusoidal) time dependent frequency, via 
the NR and coupling, was also considered. The El was 
also investigated under the same conditions. The inclu¬ 
sion of losses in the CPB and NR turns this scenario 
more realistic and the time dependence of the coupling 
X{t) and also the NR frequency uj{t) make our results 
closer to experimental conditions. The results drastically 
differ from those obtained in the resonant case. The fol¬ 
lowing scenarios were considered: (i) the resonant case 
(/ = 0)5 (^^) fhe nonresonant case, with fixed detuning, 
/ = A 7 ^ 0 , and {Hi) nonresonant case with time depen¬ 
dent detuning, f{t) = r]sm{ujft). An interesting result 
emerges: for a fixed detuning the collapse-revival effect 
does not occurs, the same being true for the El since the 
system oscillates around a value that differs from zero 
(cf. Fig. (a)) as also alerted before. However, it is 
surprising that in the case f{t) = r]sm{ujft), with the 
same conditions assumed for fixed detuning, we can see 
the El effect remaining, even in the presence of decay 
(cf. Fig. [!|(b), Fig. [7| (c)). This behavior is not shown 
with fixed detuning. In summary, we have shown that 
the use of a (time-dependent) modified BS model, which 
is extended to a more realistic scenario where the influ¬ 
ence of the losses is considered, new interesting findings 
emerge. They also indicate that it is possible to perform 
a dynamic control of the system properties by chang¬ 
ing the parameters involved. Convenient choices of the 
frequency modulation can be made to manipulate en¬ 
vironmental noisy and inaccuracies, including potential 
applications in the dynamical control of quantum infor¬ 
mation processes. We hope that these results can offer a 
reference to put the issue with force. 
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